Electronic structure, phonon spectra and electron-phonon interaction in HfB2 
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The electronic structure, Fermi surface, angle dependence of the cyclotron masses and extremal 
cross sections of the Fermi surface, phonon spectra, electron-phonon Eliashberg and transport spec- 
tral functions, temperature dependence of electrical resistivity of the HfB2 diboride were investigated 
from first principles using the fully relativistic and full potential linear mufSn-tin orbital methods. 
The calculations of the dynamic matrix were carried out within the framework of the linear response 
theory. A good agreement with experimental data of electron-phonon spectral functions, electrical 
resistivity, cyclotron masses and extremal cross sections of the Fermi surface was achieved. 
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I. INTRODUCTION 

Ceramics based on transition metal borides, nitrides, 
and carbides have extremely high melting points (>2500 
°C) and are referred to as uhra-high temperature 
ceramicsi^ Among them, diborides such as ZrB2 and 
HfB2 have a unique combination of mechanical and phys- 
ical properties: high melting points (>3000 °C); high 
thermal and electrical conductivity; chemical inertness 
against molten metals; great thermal shock resistanceii"— 
Thus, although carbides typically have the highest melt- 
ing points (>3500 °C), the diborides ZrB2 and HfB2 are 
more attractive candidates for high-temperature ther- 
momechanical structural applications at temperatures 
>3000 °Cfii^ Potential applications include thermal pro- 
tective structures for leading-edge parts on hypersonic 
re-entry space vehicles j^i^ propulsion systems jii^ fur- 
nace elements)^ refractory crucibles^ and plasma-arc 
electrodesi^ 

The discovery of superconductivity in MgB2 at 39 K 
by Akimitsu'' has lead to booming activity in the physics 
community and activated a search for superconductivity 
in other diborides. Natural candidates for this search are 
AB2-type light metal diborides (A = Li, Be, Al). How- 
ever, up to now superconductivity has not been reported 
in the majority of these compoundsi^ Only very recently 
has superconductivity below 1 K [T^ = 0.72 K) been 
reported in BeB2.75.— According to Ref. [13 no supercon- 
ducting transition down to 0.42 K has been observed in 
powders of diborides of transition metals (A = Ti, Zr, Hf, 
V, Ta, Cr, Mo, U). Only NbB2 is expected to supercon- 
duct with a rather low transition temperature (< IK), 
and contradictory reports about superconductivity up to 
rc=9.5 K in TaB2 can be found in Ref. ^ Finally, the 
reported Tc=7 K in ZrB2 encourages further studies of 
these diborides.^ 

Presently, a number of experimental studies exist deal- 
ing with the physical properties of HfB2 such as thermal 
and electrical properties r^"— mechanicalj-i^ and elas- 
tic propertiesr^ the de Haas-van Alphen (dHvA) mea- 
surements of the Fermi surface,— optical ellipsometry 
measurements magnetic susceptibilityr^i^ and NMR 



measurementSi^l First-principles calculations of the elec- 
tronic structure of diborides including HfB2 have been 
also presented] ^^'^"i^^" — 

Lawson et al^ studied the electronic structure and lat- 
tice properties of HfB2 and ZrB2 in a frame of the density 
functional theory (DFT). Lattice constants and elastic 
constants were determined. Computations of the elec- 
tronic density of states, band structure, electron localiza- 
tion function, etc. show the diverse bonding types that 
exist in these materials. They also suggest the connection 
between the electronic structure and the superior me- 
chanical properties. Lattice dynamical effects were con- 
sidered, including phonon dispersions, vibrational densi- 
ties of states, and specific heat curves. Point defect (va- 
cancies and antisites) structures and energetics are also 
presented. Vajeeston et al^ investigated the electronic 
structure of HfB2 using the tight-bonding linear muffin- 
tin orbital method, they claimed that metal-metal and 
metal-boron interactions are less significant than the p—p 
covalent interaction of boron atoms. The bonding nature, 
elastic property and hardness were investigated by Zhang 
et al}^ for HfB2 as well as ZrB2 using the first princi- 
ples total-energy plane-wave pseudopotential (PW-PP) 
method. They also reported the elastic anisotropy, Pois- 
sons ratio, hardness and Debye temperature in HfB2 and 
ZrB2. Deligoz et al.^^ investigated the structural param- 
eters (the lattice constants and bond length) and phonon 
dispersion relations in HfB2 and TaB2 compounds using 
the first-principles total energy calculations. The sec- 
ondary results on the temperature-dependent behavior of 
thermodynamical properties such as entropy, heat capac- 
ity, internal energy, and free energy were also presented. 
Zhang et al,"^^ investigated the ideal tensile and shear 
strengths of TiB2, ZrB2 and H£B2 by first-principles 
stress-strain calculations. Due to the nonlinearity of the 
stress response at large stains, the plastic anisotropy can- 
not be derived from elastic constants. Based on the rel- 
ative stiffness of boron hexagons, a bond length indica- 
tor was obtained to characterize the preference for basal 
or prismatic slip in diborides. Zhang et al?^'"^"^ inves- 
tigated theoretically the pressure dependence of elastic 
constants, bulk modulus and elastic anisotropy of HfB2. 
The pressure dependence of structural property shows 
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that the effect of pressure is httle on the structure of 
HfB2. They find high pressure greatly changes the pro- 
file of the density of states (DOS), but it hardly changes 
the DOS value at Fermi level. Meanwhile, the Mulliken 
population analyses are investigated. It was suggested 
that as the pressure increases, a number of charge trans- 
fer from Hf to B atoms. Through quasi-harmonic Debye 
model, the variations of the Debye temperature, heat ca- 
pacity and thermal expansion with pressure and temper- 
ature were obtained and discussed. 

Fedorchenko and Grechnev with coauthor o"'^^'^" mea- 
sured the temperature dependences of the magnetic sus- 
ceptibility X and its anisotropy Ax = X|| ~ Xi- for single 
crystals of transition-metal diborides MB2 (M — Sc, Ti, 
V , Zr, Hf) in the temperature interval 4.2 - 300 K. A 
transition into the superconducting state was not found 
in any of the diborides studied, right down to liquid- 
helium temperature. It was found that the anisotropy 
is weakly temperature-dependent, a nonmonotonic func- 
tion of the filling of the hybridized p — d conduction band. 
First-principles calculations of the electronic structure of 
diborides and the values of the paramagnetic contribu- 
tions spin and Van Vleck to their susceptibility show 
that the behavior of the magnetic anisotropy is due to 
the competition between Van Vleck paramagnetism and 
orbital diamagnetism of the conduction electrons. Au- 
thors of Ref. lIIi determined the thermal conductivity, 
thermal expansion. Young's modulus, flexural strength, 
and brittle-plastic deformation transition temperature 
for HfB2 as well as for HfCo.gs, HfCo.re, and IIfNo.92 
ceramics. The thermal conductivity of modified HfB2 
exceeded that of the other materials by a factor of 5 at 
room temperature and by a factor of 2.5 at 820 °C. The 
transition temperature of HfB2 was 1100 °C. Pure IIfB2 
was found to have a strength of 340 MPa in 4 point bend- 
ing, that was constant from room temperature to 1600 
°C, while a HfB2 + 10% HfC^^, had a higher room tem- 
perature bend strength of 440 MPa, but that dropped 
to 200 MPa at 1600 °C. The resuhs of the theoretical 
modeling suggest that IIfB2 should survive the high ther- 
mal stresses generated during the nozzle test primarily 
because of its superior thermal conductivity. Yang et 
alt^ used in situ spectroscopic ellipsometry to analyze 
HfB2 thin films. By modeling the film optical constants 
with a Drude-Lorentz model, the film thickness, surface 
roughness, and electrical resistivity were measured. By 
modeling the real-time data in terms of film thickness 
and surface roughness, the film nucleation and growth 
morphology were determined as a function of substrate 
type, substrate temperature, and precursor pressure. Li 
et al?~ studied the thermodynamics of the oxidation of 
HfB2 at temperatures of 1000, 1500, 2000, and 2500 K 
using volatility diagrams. They found that IIfB2 exhibits 
oxidation behavior similar to ZrB2. Zhang et alM- inves- 
tigated experimentally the thermal and electrical trans- 
port properties of various spark plasma-sintered IIfB2 
and ZrB2 based polycrystalline ceramics over the 298- 
700 K temperature range. Measurements of thermal dif- 



fusivity, electrical resistivity, and Hall coefficient were 
reported, as well as the derived properties of thermal 
conductivity, charge carrier density, and charge carrier 
mobility. Hall coefficients were negative confirming elec- 
trons as the dominant charge carrier. A Wiedemann- 
Franz analysis confirms the dominance of electronic con- 
tributions to heat transport. The thermal conductivity 
was found to decrease with increasing temperature. The 
properties of the Fermi surface of ScB2, ZrB2, and HfB2 
single crystals were studied by Phizhnikov et al}'' using 
the de Haas- van Alphen effect. The angular dependences 
of the frequencies of the dHvA oscillations in the planes 
(1010), (1120), and (0001) and the values of their effec- 
tive cyclotron masses were measured. The frequencies of 
the oscillations found lie in the interval (0.96 — 0.87) x 10^ 
T and the measured cyclotron masses m* lie in the range 
(0.09 - 0.87)too- 

Despite a lot of publications, there are still many open 
questions related to the electronic structure and phys- 
ical properties of HfB2 diboride. The most theoretical 
efforts were devoted to the lattice and mechanical prop- 
erties of HfB2. There is no theoretical expalanation of 
the Fermi surface as well as angle dependence of the cy- 
clotron masses and extremal cross sections of the Fermi 
surface, electron-phonon interaction and electrical resis- 
tivity in HfB2. The aim of this work is a complex inves- 
tigation of the electronic structure, Fermi sutface, angle 
dependence of the cyclotron masses and extremal cross 
sections of the Fermi surface, phonon spectra, electron- 
phonon Eliashberg and transport spectral functions, and 
temperature dependence of electrical resistivity of the 
HfB2 diboride. The paper is organized as follows. Sec- 
tion [H] presents the details of the calculations. Section 
mil is devoted to the electronic structure as well as the 
Fermi surface, angle dependence of the cyclotron masses 
and extremal cross sections of the Fermi surface, phonon 
spectra, electron-phonon interaction and electrical resis- 
tivity using the fully relativistic and full potenrial LMTO 
band structure methods. The results are compared with 
available experimental data. Finally, the results are sum- 
marized in Sec. IIVI 



II. COMPUTATIONAL DETAILS 

Most known transition-metal (M) diborides MB2 are 
formed by group III- VI transition elements (Sc, Ti, Zr, 
Hf, V, Nb, and others) and have a layered hexagonal C32 
structure of the AlB2-type with the space group symme- 
try PQ/mmm (number 191). It is simply a hexagonal 
lattice in which closely-packed transition metal layers are 
present alternative with graphite-like B layers (Fig. [T]). 
These diborides cannot be exactly layered compounds 
because the inter-layer interaction is strong even though 
the M layers alternate with the B layers in their crys- 
tal structure. The boron atoms lie on the corners of 
hexagons with the three nearest neighbor boron atoms 
in each plane. The M atoms lie directly in the centers of 
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FIG. 2: (Color online) Energy band structure and total DOS 
[in states/ (cell eV)] of HfB2. 



FIG. 1: (Color online) Schematic representation of the HfB2 
crystal structure. 



each boron hexagon, but midway between adjacent boron 
layers. Each transition metal atom has twelve nearest 
neighbor B atoms and eight nearest neighbor transition 
metal atoms (six are on the metal plane and two out of 
the metal plane). There is one formula unit per primi- 
tive cell and the crystal has simple hexagonal symmetry 
(D6h). By choosing appropriate primitive lattice vectors, 
the atoms are positioned at Hf (0,0,0), B (i, i, i), and B 
(|, i, i) in the unit cell. The distance between Hf-Hf is 
equal to c. This structure is quite close packed, and can 
be coped with efficiently and accurately by the atomic 
sphere approximation method. However, for precise cal- 
culation of the phonon spectra and electron-phonon in- 
teraction, a full potential approximation should be used. 

The Eliashberg function (the spectral function of the 
electron-phonon interaction) expressed in terms of the 



phonon linewidths has the for 
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27rA^(eF) Wqi. 

The line-widths characterize the partial contribution 
of each phonon: 



were N(eF-) is the electron density of states per atom 
and per spin on the Fermi level (e^) and <?k+qj'kj 
electron-phonon interaction matrix element. The double 
summation over Fermi surface in Eq.([5]) was carried out 
on dense mesh (793 point in the irreducible part of the 
BZ) 

Calculations of the electronic structure and physical 
properties of the HfB2 diborides were performed using 
fully relativistic LMTO method-"'^ with the experimen- 
tally observed lattice constants: a=3.141 A and c—iAl 
A for HfB2.'^^ For the calculation of the phonon spectra 
and electron-phonon interaction a scalar relativistic FP- 
LMTO method'^'^ was used. In our calculations we used 
the Perdew-Wang'^'' parameterization of the exchange- 
correlation potential in general gradient approximation. 
BZ integrations were performed using the improved tetra- 
hedron method. '^^ Phonon spectra and electron-phonon 
matrix elements were calculated for 50 points in the irre- 
ducible part of the BZ using the linear response scheme 
developed by Savrasov^S The 5s and hp semi-core states 
of HfB2 were treated as valence states in separate energy 
windows. Variations in charge density and potential were 
expanded in spherical harmonics inside the MT sphere 
as well as 2894 plane waves in the interstitial area with 
88.57 Ry cut-off energy for HfB2. As for the area inside 
the MT spheres, we used Zk—spd LMTO basis set energy 
(-0.1, -1, -2.5 Ry) with one-center expansions inside the 
MT-spheres performed up to Imax = 6. 



(2) 

The electron-phonon interaction constant is defined as: 

Xe-pk = 2 / —a^Fiuj). (3) 
Jo ^ 

It can also be expressed in terms of the phonons line- 
ividths: 
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7r7V(eF)wL' 



(4) 



III. RESULTS AND DISCUSSION 

A. Energy band structure 

Figure [2] presents the energy band structure and total 
density of states of HfB2. The partial DOSs HfB2 are 
shown in Fig. [31 Our results for the electronic structure 
of IIfB2 are in agreement with earlier calculations i^SiSSii^ 
A common feature for all transition metal diborides is the 
deep DOS minimum (pseudo-gap) at the Fermi energy 
separating the valence band and the conduction band. 
According to Pasturel et al.^ a pseudo-gap arises be- 
cause of a strong chemical interaction. The M-B covalent 
bonding is believed to be responsible for this effect. The 
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FIG. 3: (Color online) Partial DOSs [in states/(atom eV)] of 
Hffla. 



Hf 4f5 7 states in HIBt are situated at the —14.5 eV to 

2 ! 2 

— 10 eV. The Hf 5d states are the dominant features in 
the interval from —12.5 eV to 14 eV. These tightly bound 
states show overlap with B 2p and, to a lesser extent, with 
B 2s states both above and below ep, implying consid- 
erable covalency. Higher-energy states between 9 eV and 
17 eV above er appear to arise from Hf 6p and 6s states 
hybridized with B 2p states. The crystal field at the Hf 
site {D6h point symmetry) causes the splitting of Hf d 
orbitals into a singlet aig {d^^^-i) and two doublets eig 
{dyz and d^z) and e2g {dxy and dx2_y2). The crystal field 
at the B site {Dih point symmetry) causes the splitting 
of B p orbitals into a singlet (pz) and a doublet 62 {px 
and Py). B s states occupy a bottom of valence band be- 
tween — 14.6 eV and —3.0 eV and hybridize strongly with 
B Px and Py and Hf dyz and dxz states located at —12.5 
eV to —0.5 eV. B px and py occupied states are located 
between —12.5 eV and —0.5 eV. B pz states occupied a 
smaller energy interval from —7.5 eV to —0.5 eV with a 
very strong and narrow peak structure at around —4 eV. 



B. Fermi surface 

The Fermi surfaces (FS) of ScBz, HfBs and HfBz, were 
studied by Pluzhnikov et alM- using the dHvA effect. 

Theoretical calculations show a ring-like electron FS 
around the K symmetry point (Fig. 2]) and of a wrinkled 




FIG. 5: The calculated hole sheets of the Fermi surface at the 
A symmetry point from the 5th energy band of HfB2. 

dumbell-like hole FS at the A point (Fig. [5]) in HfBa. 
The electron and hole Fermi surfaces have threefold and 
sixfold symmetries, respectively. Figure [S] shows the cal- 
culated cross section areas in the plane perpendicular z 
direction and crossed A symmetry point for hole FS (up- 
per panel) and crossed F point for electron FS (lower 
panel) of HfB2. 
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FIG. 6: (Color online) The calculated cross sections in the 
plane perpendicular z direction and crossed A symmetry point 
(upper panel) and F point (lower panel) for HfB2 (full red 
curves) . 



Figure [7] represents angular variations of the experi- 
mentally measured dHvA frequencies^^ for HfB2 in com- 
parison with the first-principle calculations for field di- 
rection in the (1010), (1120), and (0001) planes. The 
observed frequencies of a, /3, 7, and b oscillations belong 
to electron FS around the K point. The e, /i, and C, 
orbits belong to the hole wrinkled dumbbell FS. The a 
frequencies have four branches at the (1010) plane and 
three branches at the (1120) plane. The lower 7 fre- 
quencies have one branch in both the planes. The the- 
ory reasonably well reproduces the frequencies measured 
experimentally. However, there are still some discrepan- 
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FIG. 7: (Color online) The calculated (open red and blue cir- 
cles for the electron and hole surfaces, respectively) and exper- 
imentally measuredii (black full squares) angular dependence 
of the dHvA oscillation frequencies in the compound H£B2 . 



cies. The /3 orbits have an additional two branches at 
higher frequencies at the (1120), and (0001) planes not 
observed experimentally. The experiment for high fre- 
quencies detected only e orbits in vicinity of the < 0001 > 
direction in HfB2. We found the e, /i and <^ orbits simi- 
lar to the corresponding orbits observed experimentally 
in isostructural and isovalent ZrB2.^^ These orbits have 
not been detected in the dHvA experiments^ One of the 
possible reasons for that is the relatively large cyclotron 
masses for these orbits. Figure [8] shows the theoretically 
calculated angular dependence of the cyclotron masses 
(jjih) and the experimentally measured masses (m*) for 
high symmetry directions in HfB2. The cyclotron effec- 
tive masses were determined from the temperature de- 
pendences of the amplitudes of the dHvA oscillations. 
The cyclotron masses for the e, /i, and C orbits in HfB2 
are much higher than the corresponding low-frequency 
oscillations a, /3, 7 and 8. The fact that the masses for 
electron Fermi surface are significantly larger than for the 
hole Fermi surface may explain a negative experimentally 
measured Hall coefficient^i^ and confirms electrons as the 
dominant charge carriers in HfB2. A Wiedemann- Franz 
analysis also indicate the dominance of electronic contri- 
butions to heat transport 

We note that band cyclotron effective masses mf, arc 
renormalized by the electron-phonon interaction m* = 
mb(l -|- A), where A is the constant of the electron-phonon 
interaction. By comparing the experimentally measured 
cyclotron masses with band masses we can estimate the 
A. It is strongly varied on the orbit type and magnetic di- 
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FIG. 9: (Color online) Theoretically calculated phonon den- 
sity of states (full blue line) for HfB2. The dotted green and 
red dashed lines present the calculated phonon DOS of HfB2 
by Deligoz et alr^ and Lawson et alM; respectively. 



FIG. 8: (Color online) The calculated angular dependence 
of the cyclotron masses for the electron Fermi surface (open 
red circles) and the hole Fermi surface (blue open triangles) 
and experimentally measured ones^^ (black full squares) in 
the compound IIfB2. 



rection. We estimate the constant of the electron-phonon 
interaction to be equal to 0.18-0.23 for the a orbits and 
0.36 and 0.75 for the e and orbits, respectively, with 
H||< 0001 >. For the < lOlO > and < 1120 > directions 
the A for the a orbits are reduced, respectively, to 0.10 
and 0.12 values. 



C. Phonon spectra 

The unit cell of HfB2 contains three atoms, which gives 
in general case a nine phonon branches. Figure |9] shows 
theoretically calculated phonon density of state for HfB2 
(full blue curve). The DOS for HfB2 can be separated 
into three distinct regions. Based on our analysis of rela- 
tive directions of eigenvectors for each atom in unit cell, 
we find that the first region (with a peak in phonon DOS 
at 5.2 THz) is dominated by the motion of Hf. This re- 
gion belongs to the acoustic phonon modes. The second 
wide region (14-20 THz) results from the coupled motion 
of Hf and the two B atoms in the unit cell. The 
^2g, Big phonon modes (sec Table|T| lie in this area. The 
phonon DOS in the third region extends from 22 THz to 
26 THz. This is due to the movement of boron atoms and 
is expected since boron is lighter than Hf. The covalent 
character of the B-B bonding is also crucial for the high 
frequency of phonons. The in-plane Ei2g mode belongs 
to this region. The second and third regions represent 
optical phonon modes in crystals. The most significant 



TABLE I: Theoretically calculated phonon frequencies (in 
THz) in the F symmetry point for HfB2 and calculated 
phonon frequencies by Deligoz et al^ and Lawson et al. Ref. 



reference 


E\u 




Big 


E2g 


our results 


13.76 


15.03 


17.12 


25.17 


SIESTA2S, 


14.10 


15.19 


15.87 


24.49 


VASP^a 


13.34 


14.00 


16.40 


24.16 


ABINITi22 


12.92 


13.85 


16.01 


23.59 



feature in the phonon DOS is a gap around 6-13 ThZ. 
This gap is a consequence of the large mass difference 
between B(10.8 a.u.) and Hf (178.49 a.u.), which leads 
to decoupling of the transition metal and boron vibra- 
tions. 

Currently, there are no data concerning the experimen- 
tally measured phonon DOS in HfB2. So we compare 
our results with theoretically calculated phonon DOS by 
Deligoz et al.^^ and Lawson et al^ (Fig. M and Ta- 
ble HI). Calculations of Deligoz et al^ were based on 
the density functional formalism and generalized gradient 
approximation. They used the Perdew-Burke-Ernzerhof 
functional-? for the exchange-correlation energy as it is 
implemented in the SIESTA codeM- This code calcu- 
lates the total energies and atomic Hellmann-Feynman 
forces using a linear combination of atomic orbitals as 
the basis set. The basis set consists of finite range pseu- 
doatomic orbitals of the Sankey-Niklewsky type**- gener- 
alized to include multiplezeta decays. The interactions 
between electrons and core ions are simulated with the 
separable TrouUier-Martins^ normconserving pseudopo- 
tentials. In other words, they used the so-called "frozen 
phonon" technique and built an optimized rhombohedral 
supercell with 36 atoms. This method is inconvenient 
for calculating phonon spectra for small q-points as well 
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as for compounds with large number of atoms per unit 
cell. Lawson et al^ used two different codes to calcu- 
lated the phonon spectra. VASP, the supercell method, 
based un the projected augmented wave potentials. Sec- 
ond method, ABINIT, used Fritz Haber Institute pseu- 
dopotentials in the Troulliers-Martin form. VASP results 
of Lawson et al.^^ is slightly closer to our calculation with 
respect to ABINIT data. There is a good agreement be- 
tween our calculations and the results of Deligoz et al^ 
in a shape and energy position of the second peak in the 
phonon DOS. There is an energy shift towards smaller 
energies of the first and third peaks of the Lawson et 
al^ calculations in comparison with the Deligoz et 
data with our results are just in between these two cal 
culations. 
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D. Electron-phonon interaction 

Figure [TU] shows theoretically calculated Eliashberg 
functions for IIfB2 as well as electron-phonon prefac- 
tor a^(w) (definition of this function is merely ratio 

(uj) F (uj) / F (uj)) . There is no difference between main 
peaks positions of phonon spectra and electron-phonon 
coupling function. Electron-phonon prefactor has three 
peaks: 5.2 THz, 17.1 THz and 21.3 THz (the corre- 
sponding peaks in the phonon DOS are situated at the 
5.2 THz, 17.1 THz, and 23.4 THz frequencies). The 
a^(a;) has strongly varying character. Therefore the 
electron-phonon coupling can not be factorized into in- 
dependent electronic and phonon parts. The matrix el- 
ement of electron-phonon interaction cannot be repre- 
sented in form a^{Lu) « const and hence well known 
McMillan approximation"*^ is not valid for HfB2. By inte- 
grating the Eliashberg function we estimate the averaged 
electron-phonon constants Ae_pft=0.17. The constant of 
the electron-phonon interaction also can be estimated 
by comparison the theoretically calculated DOS at the 
Fermi level with the electron specific heat coefficient 7. 
Cp = -fT, where 7 =1.0 mJmole-^K-^ for HfBa^i HfBa 
possesses quite small value of the DOS at the Fermi level 
of 0.4 states/(cell eV), it gives the theoretically calcu- 
lated 7b=0.8 mJmole~^K~^ and A=0.2 with qualitative 
agreement with Ae-p/i=0.17. 



FIG. 10: (Color online) Theoretically calculated Eliashberg 
function a^F{uj) of HfB2 (full red line) and electron-phonon 
prefactor a^{uj) (dashed blue line). 



where, the subscript / specifies the direction of the elec- 
trical current. In our work, we investigate two direc- 
tion: [0001] (c-axis or z direction) and [1010] (a-axis or 
x-direction) . (w|) is the average square of the / compo- 
nent of the Fermi velocity, ^ = uj/2kBT. 

Mathematically, the transport function atrF{uj) dif- 
fers from aFiuj) only by an additional factor [1 — 
vi{'k)vi{'k')/ (vj)], which preferentially weights the 
backscattering processes. 

Formula ([5]) remains valid in the range Ot, /5 < T < 
2Qtr^^ where: 



2 



,2\l/2 



ujaf^F{uj)du}, 



Xtr — 2 



alFi^) — , 



, duj 

UJ 



(6) 



(7) 



(8) 



E. Electrical resistivity 

In the pure metals (excluding low-temperature region), 
the electron-phonon interaction is the dominant factor 
governing electrical conductivity of the substance. Using 
lowest-order variational approximation, the solution for 
the Boltzmann equation gives the following formula for 
the temperature dependence of pi{T): 



Pi{T) = 



TTV,cellkBT [°° duj ^ 



rF{L0), (5) 



The low-temperature electrical resistivity is the result 
of electron-electron interaction, size effects, scattering on 
impurities, etc., however, for high temperatures it is nec- 
essarily to take into account the effects of anharmonicity 
and the temperature smearing of the Fermi surface. In 
our calculations Otr=654.4 K for c-axis, and 679.9 for 
a-axis for HfB2. 

Figure [TT] represents the theoretically calculated tem- 
perature dependence of electrical resistivity of HfB2 for 
the < 0001 > direction (full blue curve) and the basal 
< 1010 > direction (dashed red curve) and experimen- 
tal measurements for polycrystalline HfB2.— Specimen 
of ceramic HB2 was obtained by spark plasma sintering 
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FIG. 11: (Color online) Theoretically calculated for the < 
0001 > direction (full blue curve) and the basal < lOlO > di- 
rection (dashed red curve) and experimentally measured tem- 
perature dependence of electrical resistivity of HfB2.*'* 

method and had good ratio of experimental and theo- 
retically calculated density Pexplpth = 98.1% Our re- 
sults are in good agreement with the experiment. The 
small discrepancy does not exceed accuracy of calcula- 
tion. We obtained anisotropy ratio of electrical resistivity 
at T=300K: Pz/ Px = 1.079. Actually this fact indicates 
that for HfB2 anisotropy is not clearly expressed. 

IV. SUMMARY 

We have studied the electronic structure and physical 
properties of HfB2 using fully relativistic and full po- 
tential linear muffin-tin orbital methods. We study the 
electron and phonon subsystems as well as the electron- 
phonon interaction in this compound. 

We investigated the Fermi surface, angle dependence 
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of the cyclotron masses, and extremal cross sections of 
the Fermi surface of HfB2 in details. Theoretical calcula- 
tions show a ring-like electron FS in HfB2 around the K 
symmetry point and a wrinkled dumbbell-like hole FS at 
the A point. Theory reproduces the experimentally mea- 
sured dHvA frequencies in HfB2 reasonably well. We 
found that masses for low-frequency oscillations a, (5, 7, 
and 5 are less than 0.2577io. Masses for high-frequency os- 
cillations e, p., and C he in the range from —0.3 to —0.65 
toq. The experiment for high frequencies detected only e 
orbits in vicinity of the < 0001 > direction in HfB2. We 
found the e, pL and C, orbits similar to the corresponding 
orbits observed experimentally in isostructural and iso- 
valent ZrB2. These orbits have not been detected in the 
dHvA experiment. One of the possible reasons for that 
is the relatively large cyclotron masses for these orbits. 

Calculated phonon spectra and phonon DOSs for HfB2 
is in good agreement with previous calculations. We 
did not found regions with high electron-phonon inter- 
action or phonon dispersion curves with soft modes in 
HfB2. This is in agreement with the fact that no trace 
of superconductivity was found in these borides. The av- 
eraged electron-phonon interaction constant was found 
to be rather small Ae-p/i=0.17 for HfB2. We calculated 
the temperature dependence of the electrical resistivity in 
HfB2 in the lowest-order variational approximation of the 
Boltzmann equation. We found rather small anisotropic 
behavior of the electrical resistivity in HfB2 to be in good 
agreement with experimental observation. 
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